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INTRODUCTION 


2 


The quadratic equation zf — x — 1 = 0 is known to have roots 


EVS 
ee 
but suppose that we decide instead to solve this equation using an iterative 
method. The equation can be rearranged as 


z=1+-, 
£ 


so we could find successive approximations, £n, to a root by using the 
formula 


1 

In4+1 = 14+—. 

In 

If we take x; = 1 as our initial approximation, the sequence begins as 


follows. 


Dia 
1 2 
T2 eS +7 I 
1 1 3 
r3=1+—=1+ =5 
2 Tyee 
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It is interesting to see the Fibonacci numbers appearing as numerators and 
F, n+1 


denominators of each £n, and we might conjecture that £n = , for all 


n 
positive integers n. This is certainly true when n = 1, and if we assume that 


F, 
te= ra then 
1 F; Fri + F; F; 
P E E ee 
Tk Fra Fk+1 Fk+1 
So Mathematical Induction confirms that it is true that £n = N , for all 
n 


Mo. 


The question remains as to whether or not the sequence 21, £2, £3, ... 
produced by this iterative process converges. The question is the same as 
Fo F} F, 

Bi Po aie 
Fibonacci numbers converges, and this problem has a nice geometric 
interpretation. Imagine squares with edge lengths 1, 1, 2, 3, 5, 8, 13, .... 
These can be joined, in order, to form a sequence of rectangles where, at 
each stage, the new rectangle is formed by adjoining the next square to the 
longer edge of the existing rectangle. This process is illustrated in the 
following figure. Now the shape of each rectangle, as determined by the ratio 
of the longer to the shorter side, is the ratio of successive Fibonacci 
numbers. So the question becomes: does this sequence of rectangles have a 
limiting shape? 


asking whether the sequence of ratios of successive 


The roots of ax” + bx + c = 0 are 


gf —b + Vb? — 4ac 
a 


We met the Fibonacci numbers in 
Section 5 of Unit 2. 


hex]: 


5x5 
3x3 


13 x 13 


8x8 


The answer, as we shall see from the work of this unit, is yes. If, for the 
moment, we assume that the sequence has a limit k, say, then it is not too 
difficult to see what k must be. For sufficiently large rectangles the ratio of 
the lengths of edges is very close to k, so we could take these lengths to be 
(approximately) ky and y. To obtain the next rectangle we add on a square 
of side ky as shown below. 


ky ky 


The ratio of the sides of this larger rectangle is also (approximately) k and so 


ytky _ k, 
ky 
which, after cancelling the y, simplifies to k? — k — 1 = 0. Hence, if the 
sequence tends to a limit, it will be a root of the quadratic equation from 
l+v5 | 


2 
number which is known as the golden ratio. The value of the golden ratio is 


approximately 1.6180. 


which we started, z? — x — 1 = 0. The positive root is k = 


The golden ratio was, along with m and v2, one of the first irrational 
numbers known to, and used by, the ancient Greeks, arising originally in the 
geometrical problem of dividing a line segment AB at a point C such that 


bas = this common ratio turning out to be the golden rati 
AC ~ GB’ s o io turning golde io. 
A C B 
-—— — — 


The golden ratio crops up all over mathematics and is particularly 
prominent in problems concerning regular plane and solid figures. 


As one example, it is the ratio of the radius of a circle to the edge-length of 


an inscribed regular decagon. The ancient Greeks appear to have been Trigonometry shows that this ratio 
endeared with the golden ratio and it may have been deployed by their is 1 
sculptors and architects for aesthetic purposes. It features in the structure of 2sin(7/10)’ 


the Parthenon and the theatre at Epidavros; for example, the theatre has 21 which is another expression for the 
tiers of seats above the gangway and 34 below, 21 and 34 being consecutive golden ratio. 
Fibonacci numbers whose ratio is a good approximation to the golden ratio. 


Our main purpose in giving this example is to illustrate what is meant by a 
continued fraction. Each of the expressions 


i ea ee ee ae = 
t+ I 1+ ER DES o 
1+ I 1+ tT 
I+ i 
is called a finite continued fraction. Their limit is the infinite continued 
fraction: 
I : I 
14+ 
1 
L4 I 


This is the simplest, and best-known example of a continued fraction. In the 
more general infinite continued fraction 


a, + 
a2 + 


e a 


the a; are positive integers, with the exception of the first one, a1, which 
may be a positive or negative integer or zero. 


The theory of continued fractions is extensive and their area of applicability 
very wide, so in this unit we can only hope to give you something of the 
flavour of them. We shall give a number of diverse applications and, in 

Unit 8, we shall see how they can be exploited to lead to a neat solution of a 
famous Diophantine equation. 


We shall first explore finite continued fractions and, in so doing, lay the 
foundations for our subsequent investigations of infinite continued fractions. 
We shall discover that every infinite continued fraction converges to an 
irrational number and, conversely, that every irrational number can be 
expressed, in a unique way, as an infinite continued fraction. 


We are all in the habit of writing down, and even manipulating, expressions 


14+ V5 


when we come to practical calculations involving them we must deal with 
rational numbers. (Even today’s computers are limited in how many decimal 
places of 7 they can handle.) So the question of finding good rational 
approximations to irrational numbers is an important one. In the final 
section of the unit, we see how the continued fraction of an irrational 
number provides us with ‘best possible’ rational approximations. 


such as 7, log, 3, , etc., which represent irrational numbers. But 


1 FINITE CONTINUED FRACTIONS 


1.1 The continued fraction of a rational number 
We start with the key definitions. 


Definition 1.1 Finite Continued Fractions 


Let a; be any real number and let a2, a3,...,@n, be any positive real 
numbers. The expression 


1 
a + 


a2 + 
a3 + 


Ot DDE 
an 


is a finite continued fraction. The numbers a1, a2, ..., an are the 
partial quotients of the continued fraction. When all the partial 
quotients are integers the continued fraction is said to be simple. 


Continued fractions can be cumbersome objects and it is therefore 
advantageous to have a convenient symbol to display them. Any continued 
fraction is determined by its partial quotients and so we write the finite 
continued fraction in the definition as 


[a1, a2, 03;... 1 An]. 


For example, 
1 
9 pa, 
T 5 
; y aotea ; 1G es 
and the expression on the right simplifies to the rational number li which is 


the value of the continued fraction. Similarly the value of [—3, 1,2, 4] can be 
calculated as follows. 


[—3, 1, 2,4] = —3 + nea =—3 + = = —3 + £ = a Note that formally we think of the 
= Jes 13 13 finite continued fraction and its 
24 1 9 value as being distinct. The 
4 continued fraction is an arithmetic 
representation of a number (its 
Problem 1.1 value). However we shall 
? : , P occasionally drop this distinction 
Determine the rational numbers represented by the following finite and write, for example, 
continued fractions. x = [—3, 1,2, 4] to mean that x is 
(a) [2,3,4] (b) [0,2,2,2] the value of this continued fraction. 


Looking at the first step in the calculations in the examples above, note that 


1 1 
1,2,5] = 1 + — d [-3,1,2,4] = —3 + —. 
a =) aang poe eal 
In general, from the way the finite continued fraction is defined, we have the 
following property. 


First continued fraction identity 
1 


[a2,3,.--, An] 


[a1, 2, @3,.--,An| =a,+ 


In our work we are going to be concerned almost exclusively with simple 
continued fractions; that is, ones in which all the partial quotients are 
integers. There will, however, be isolated occasions when it proves 
convenient to use continued fractions with real partial quotients for the 
purpose of providing general proofs, though the results will subsequently 
only be applied to cases involving simple continued fractions. 


In the text, for ease of notation, we shall use the abbreviation FCF to stand 
for a simple finite continued fraction. 


It is readily seen from the definition that the value of any FCF is a rational 
number. As it happens the converse also holds: any rational number can be 
expressed as an FCF. The Euclidean Algorithm provides a means of 
determining an FCF for a given rational number. We illustrate this method 
in the following example. 


Example 1.1 
225 


D i F —. 
etermine an FCF for 157 


Dividing 225 by 157, the Euclidean Algorithm gives the sequence of 
equations on the left below. Each of these is written in an alternative form 
on the right. 


295 = 1x 157+68 or = a — 
157 = 2x 68421 D 
68 = 3K 5 S ja 
91 = 4x54+1 = 44s 

5 = 5x1 *=5 


From the system of equations on the right, by continually using the next 
equation to substitute for the existing fraction, we obtain the following. 


ce r 1 
157 Lak of 157 
68 
1 1 
=1+ z ~ 1+ T 
T 2+ 55 
21 
1 1 
=1+ =1+ 
1 1 
2+ 5 2+ T 
3+ 57 ‘iar 
5 
=1+ i 
= 1 
2+ I 
eae 
4-- = 
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225 


In other words, is7 [1, 2,3,4,5], the partial quotients of the FCF being, in 
order, precisely the quotients as they arise in the Euclidean Algorithm. + 
Problem 1.2 

Determine an FCF for each of the following rational numbers. 


aF ws 


This way of using the Euclidean Algorithm can be deployed to obtain an 
FCF for any rational number. Let us formalize this before trying more 
examples. 


Theorem 1.1 Equivalence of FCFs and rational numbers 


Every simple finite continued fraction has a rational value and, 
conversely, every rational number can be expressed as a simple finite 
continued fraction. 


Proof of Theorem 1.1 


We prove that the value of any FCF is rational by Mathematical Induction 
on the number of partial quotients. 


The induction is started by observing that any FCF with a single partial 
quotient, [a1], is an integer and therefore rational. By definition [a1] = a1. 


For the induction step, suppose that the value of any FCF with k partial 
quotients is rational and consider the FCF [a1, a2, a3,...,@%41| with k +1 
partial quotients. We now use the first continued fraction identity: 


1 


Q1, 42,03,...,QAk41] = @y + a. 
| r ; — | [a2, a3, . . - , ak+1] 

Now a1, being an integer, is rational. Furthermore, by the induction 

hypothesis, [a2, a3, ...,@41] is rational and hence so is its reciprocal. It As az > 0, by definition, 
follows that [a1,a@2,a3,...,@%41], being the sum of two rationals, is itself [az,@3,.-.,@%41] £0. 
rational. 


Mathematical Induction therefore confirms that the value of any FCF is a 
rational number. 


a 
Conversely, let 5 be any rational number. By the Euclidean Algorithm: When we write © for a rational 
a z number we assume that a, b are 
a=Mmb+n, O<rn<b or == qF es integers with b positive and 
b b gcd(a,b) = 1. 
b T2 
b=qritre, O<T2<T1 or — =@+— 
ry ry 
Ti T3 
ry = qQ3r2 +r3, O<T3<T2 or —=q+— 
T2 r2 
Tn- Tn-1 
Tn-3 = In-17n-2 +Tn-1, O< Ta- < Tn- of == cgp H 
Tn—2 Tn—2 
= +0 or La 
Tn—2 = QnTn—1 ’ Tai an The fact that the non-zero 


remainders decrease guarantees 
that the algorithm stops after a 
finite number of steps. 


From the right-hand system of equations we obtain 


a 1 1 il 
gre | it 
a tT q2 + I 
T2 F Tz 
T3 
F 1 
at a= I , 
qit I 
q3 + 
Á 1 
neta 
dn 
which shows that S has the FCF [q1, 92, --, qn]. a 
Problem 1.3 
Obtain FCFs for the following rational numbers. 
65 2a —94 
(a) z3 b 33 O 


The FCF of a rational number is not unique. In fact each rational number 
has exactly two FCFs. The ambiguity arises because of the freedom in 
choosing whether or not to terminate the FCF with partial quotient 1. To 


substantiate that remark, look again at Example 1.1 where we saw that 
225 i 
ir [1,2,3,4,5]. The final partial quotient, 5, can be written as 4 + T 


1 ih 
piek rene 


I : 
4+ 5 4+ T 
eS 
i 1 
Consequently we have an alternative form for the FCF 


225 
= = 4,4, 1}. 
ee = 11,234,451] 


Similarly, from Problem 1.3 part (a) above 


65 
z = [21,4,1,3]= ft 44,9. i, 


1 
the alternative being reached by writing the final partial quotient 3 as 2 + rt 
In general, if an > 1, 


|a1, a2,. oe Gere = (a1, a2,. ++)Qn—-1,4n — k Ll; 


1 
since an = (an — 1) + T 


Problem 1.4 


Write each of the following as an FCF with an even number of partial 
quotients. 


(a) [1,2,3] (b) [-4,3, 1,6,1] (c) [0] 


If an FCF has a final partial quotient of 1 then, in the way illustrated in 
Problem 1.4 part (b), the final 1 can be absorbed into the previous partial 
quotient. If we choose always to absorb such 1s then the ambiguity in the 
FCF is removed and every rational number has a unique FCF whose final 
quotient is other than 1. (There is one exceptional number. The integer 1 
has FCF [1], where the final partial quotient 1 cannot be absorbed into the 
previous one because there is no previous one. The two FCFs of 1 are [1] and 


10 


(0, 1], both of which end in 1.) We shall not give a proof of this uniqueness 
here for, in Section 3, we shall prove a general result concerning the 
uniqueness of infinite continued fractions which can be applied to the finite 
case. 


1.2 The convergents of a continued fraction 


Let us turn to the question of evaluating FCFs. First we need to introduce 
some further terminology. 


Definition 1.2 Convergents 


For each k = 1,2,...,n we define the kth convergent, Cp, of the FCF Notice that the final convergent Cp 
[a1,@2,...,@n] to be the value of [a;,a2,... , Ox]. is the value of the FCF itself. 


The reason for the name convergent 
will become apparent in due course. 


As an illustration of the definition consider [1, 2,4, 3,2]. It has five 


convergents. 
C= |) =1 
Lea 
Co] T a = 4 
1 13 
Cea 52,4) = 1+ — = 
2+- . 
4 
1 42 
Cie Aes eet —— Se Although, for example, we have, 
ee ag written C4 = [1,2, 4,3] the 
A+ 1 convergent C4 is really the value of 
i 97 this expression, namely 29° 
Cx —|1,254°3> 2) a4 ae acs 
as S 
4+ ETSE 
gar 
= 2 


At first sight it appears that the only way of evaluating an FCF is the 
logical one of starting with the final partial quotient and working our way 
upward through the expression, determining the various fractions along the 
way. Fortunately this is not so. We can evaluate the convergents 

C1, C2, ... directly, obtaining each one from the convergents which have 
gone before it until we reach the one we want. Our next result will show 
how this is done. 


A key step in the proof of the coming theorem is going to be the following 
identity. 


Second continued fraction identity 


1 
[a1,2,..., A, 0x41] = [a1,Q2,...,@4—-1, aK + ] 
Qk+1 


This is obtained by absorbing the final partial quotient of the left-hand side 
into the penultimate one. This enables us to write a continued fraction with 
k + 1 partial quotients in an alternative form with just k partial quotients. 
However, the expression on the right-hand side of this equality is not a 


simple continued fraction because, in general, ag + is not an integer. ak + 


Qk41 Qk+1 
For this reason we shall state, and prove, our theorem for the case of the when ak+ı = 1. 


will be an integer only 
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general continued fraction, although thereafter our applications will be 
restricted to the special case where the continued fraction is simple. 


Theorem 1.2 The value of a finite continued fraction 


Let [a1,@2,...,@n| be any finite continued fraction. Define the 
numerators pı, P2, P3, ---, Pn and the denominators qi, q2, 43, --- 
of this continued fraction by: 


Pı = 41, p2 = a12 + 1; Pk = GkPK—1 + Pk-2, 3<k<n; 
nj =1, q2 = a2; Gk = @kQk—1 + qk-2, askin: 


Then the value of [a1,a2,...,@n| is p 


Before we tackle the proof, let us get a feel for the formulae involved here by 
working through some examples. 


Example 1.2 
Use the formulae of Theorem 1.2 to determine the value of [3,3, 1, 2, 5]. 
The first two numerators are 
pi =a; =3 and po=aja2+1=3x3+4+1=10. 
Subsequent numerators are given by the iterative formula 
Pk = AkPk-1 + Pk-2- 
p3 = 1 x p2 + pı = 1 x 10 + 3 = 13 
pa = 2 x p3 + p2 = 2 x 13 + 10 = 36 
Ds = 5 X pat p3 = 5 X 36 + 13 = 193 
The first two denominators are 
g= 1 andaga a> Sa; 
Subsequent denominators are given by the same iterative formula as for the 
numerators, qk = @kQk—1 + qk-2- 
gg =1xXgq@tqga=1x34+1=4 
q4 =2xgtq=2x44+3=11 
qs =5Xqat+q3=5x11+4=59 
193 
59° 


The value of [3,3, 1, 2,5] is, according to the theorem, 2 = + 
5 


Notice that the calculations in the above example give us more than just the 
value of the FCF; we have determined the values of all the convergents along 
the way. For example, the theorem tells us that the third convergent [3, 3, 1] 


13 
has value ees —. Indeed the convergents are precisely eh E2. P3.. pa and 
q3 š w Tye =% qı 9 93 Q4 
as That is, ae 8 ag and 9” the last of which is the value of the 


FCF itself. 


This example illustrates an immediate corollary to the above theorem. 


Corollary to Theorem 1.2 


Let [a1,@2,...,@n] be any finite continued fraction and the numerators 
P1, P2,- - -Pn and denominators qi, q2,- --,qn be as defined in the 
theorem. Then its convergents are 


tonk = Aaea 
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Problem 1.5 


Use the corollary to Theorem 1.2 to determine the convergents of the 
following FCFs. 


(a) [2,3, 4,1, 2] (b) [—2,3, 1,6] 


In the next example we suggest a tabular method for calculating the 
successive numerators and denominators. This method may help you avoid 
arithmetic slips and might well assist you in remembering these important 
formulae. However, after this subsection we shall not display tables again; 
we shall just use the formulae to write down the convergents. 


Example 1.3 
Determine the rational number represented by [4, 2, 1, 3, 1, 2, 4]. 


In the table we set out four rows labelled k, pk, az and qp. Since this 
continued fraction has seven partial quotients, k takes the values from 

1 to 7. The a, are the given partial quotients and so the entries for this row 
can be displayed immediately. The p and qj, rows are started by entering 
the first two values for each using the formulae pı = a1, p2 = aja + 1, 

qı = 1 and Q = a2. 


To complete the table we move along the a, row filling in entries above and 
below each ax. For the entry above aj, (i.e. px), we multiply a, by the entry 
diagonally above and to the left, and add to this product the entry to the 
left of the latter (Figure 1.1). So, for example, p = 1 x 9 + 4 = 13 and 

p4 = 3 x 13 +9 = 48. 


Pk = GkPr—1 + Pk-2 Ok = @kQk—1 + Qk-2 


Figure 1.1 Figure 1.2 


The same procedure applies for the q,. We multiply ap by the entry 
diagonally below and to the left and add on to this product the entry to the 
left of the latter (Figure 1.2). So, q3 = 1 x 2 + 1 = 3, and so on. 


13 


The complete table is 


4 9 13 48 61 170 741 


and the convergents are 12’? 3°’ li’ i’ 39 and 170’ the final one being 
the value of the FCF. 4 
Problem 1.6 


Use the tabular method to determine the convergents of each of the 
following FCFs. 


(a) [—4, 1, 5, 1, 6] (b) (1,1, 2,1, 2, 1,2] 


Problem 1.7 


In Problem 1.5 part (a) we determined the convergents of [2,3, 4, 1,2]. Now 
determine the convergents of the following FCFs. 


(a) [2,3,4,1,2,2]  (b) [2,3,4,1,2, 100] 


Problem 1.7 highlights another property of continued fractions. The 

cae and = and any continued 

ea 5 7 

fraction which begins [2, 3, 4,1,2,...] will have these five as its first, second, 

third, fourth and fifth convergents no matter what follows. There is nothing 

special about [2, 3,4, 1,2] here. If [ai,a2,...,@n] is any continued fraction 

with convergents = a Aer = then the continued fraction obtained by 
qi 42 n 

adding a further partial quotient, [a1, a2, .- - , an, @n+1], has the same 


convergents = p2 


convergents of [2, 3, 4, 1, 2] are 


s Z but with one extra final one, which will be equal 
1 %2 n 

to the value of the (new) continued fraction itself. According to 

Qn+1Pn + Pn-1 


Q@n419n + In-1 
have to establish for a proof by induction. 


Theorem 1.2, this final convergent is , and this is what we 


With this observation we are now ready to prove Theorem 1.2. 


Proof of Theorem 1.2 


We prove the result by Mathematical Induction on the number n of partial 
quotients in the continued fraction. 


As the basis for the induction, we note that the result holds for any 
continued fraction [a1] with a single partial quotient and for any continued 
fraction [a1, a2] with two partial quotients since 

ay 


Pı 
== d = = — 
lai] =a, an = i ay 


and 
pao aa2+1 _ 1 


@ı + —. 


1 
[a1,a@2]=a,+— and 
a2 q2 a2 a2 
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Now suppose that the result holds for any continued fraction with k > 2 
partial quotients and consider the finite continued fraction 


[a1,@2,..., ak+ı] 


with k + 1 partial quotients. The numerators of this continued fraction are 


Pi, P2, ---, Pk and ppk+1, and the denominators are qi, q2, ---, qk and qk+1- 


Our task is to show that the value of the continued fraction is Pet) 


dk+1 
Now, by the second continued fraction identity, 


1 / 
[@1,@2,...,@%41] = [a1,@2,...,@x%—1,a%], where aj, = ag + i 
k+1 


Look at the right-hand continued fraction. Its numerators are 


P1,P2,---,Pk-1 and p, = (a F ) pes + Pr-2- 


Qk4+1 


Similarly its denominators are 


1 
91;92,-+-;9k-1 and gq, = (a zy >) dk-1 + qk-2. 
Ak+1 


Moreover, as it has just k partial quotients, the induction hypothesis tells us 
/ 


that its value is Pr. That is, 
dk 


pi, 
lepar- sagi] = ai 
dk 


f 
(a ap ) Pk—1 + Pr—2 
Ak+1 


1 
(a a ) Qk—-1 + Qk-2 
Qk+1 


_ Oe 41(aKPR—1 + Pe—2) + PR-1 
ak+ı(akqk—1 + qk-2) + qk-1 

_ ûk+1Pk + Pk-1 _ Pk+1 

 k+19k + qk-1 k+” 


This completes the induction step and therefore, by Mathematical 
Induction, the result is true for a continued fraction with any number of 
partial quotients. a 


1.3 Properties of the convergents 


The convergents of an FCF have a number of curious properties. For 
example, let us look again at the sequence of convergents of [4, 2, 1,3, 1, 2, 4] 
which we found in Example 1.3. 


m 4 9 13 48 61 17 741 
oto” Sy a 39 "170 


The first thing you may have noticed is that, for each k, the numbers pọ and 
qk are relatively prime. There is also something curious about the 


‘ ? k- 3 
cross-products’ of consecutive convergents. For example, for — and — we 


have 9 x 3 = 27 and 13 x 2 = 26. They differ by 1, something which 
happens with every pair of consecutive convergents in this example. You 
might like to check that 170 x 14 — 61 x 39 = 1. This is a general property 
of the convergents of an FCF. 
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Theorem 1.3 Properties of the convergents 


Let [a1, a2, ...,an] be any FCF, and let pı, po, ..., Pn and 
q1, G2, ---, Mn be its numerators and denominators as defined in 
Theorem 1.2. Then we have the following properties. 
(a) Pkdk—1 — Pk-14k = (—1)*; that is, 
Z —1)* 

Pe Pt for2<k<n. 

dk = Qk-1 GkWk-1 
(b) prak—2 — pk-2qk = (—1)*~1a,; that is, 


Theorem 1.2 told us that the 
convergents are Te; Property (c) 


-2 CD 
Pe a SN ee E 
dk qk—-2 qk Vk—2 


confirms that these numbers are 
rationals, written in their lowest 


(c) pk and gx are relatively prime for k = 1, 2,..., n. oe, 


Proof of Theorem 1.3 


(a) We use Mathematical Induction to prove that, for any integer r > 2 the 
numerators pr—ı and p, and the denominators qr—ı and qr of the FCF 


[a1,@2,...,a,]| satisfy 
PrOr—1 — Pr—1gr = (-1)’. The result of property (a) will 
follow by letting r range from 
For the case r = 2, pı = a1, po = a1đa2 +1, qı = 1 and q2 = az so that 2ton 


pod — Pig2 = (aaz + 1) x 1 — a102 = 1 = (-1)?. 
This gives the basis for the induction. 
For the induction step, assume that it is true for r = k, so that, for the 
FCF [a1, a2,...,@%], we have that prqx—1 — pk-1qk = (—1)*, and 
consider |[a1, a2, ..., ak, @ķ+1]. Its final numerator and denominator are, 
according to Theorem 1.2, 
Pk+1 = @k4+1Pk+Pr-1 and qk+1 = e419 + qk-1- 
Therefore 
Pk+i9k — Pkak+1 = (Ak+1Pk + Pk-1)qk — Pk(Ak419k + gk-1) 
= Pk-19k — PkQk-1 
= (—1)(Pkqk-1 — Pr—19k) 
=(-—1)(—1)*, by the induction hypothesis, 
= (1)! 
showing property (a) to be true when r = k + 1. This completes the 
proof. 


(b) This is a simple consequence of property (a). 


PkOk—2 — Pk—29k = (AkPk—1 + Pk—2)k—2 — Pk-2(akqk-1 + qk-2) 
= Ak (Pk—19k—2 — Pk-24k—1) 
= ax(—1)*"! 
(c) This also follows from property (a). 


If d = gced(px, qk) then prqx—1 — Pk—19k, being an integer combination of 
pr and qk, is a multiple of d. Therefore, by property (a), d divides 
(—1)*, which implies d = 1. a 


Let us pause at this point to look at one application of FCFs. The result of 
Theorem 1.3 leads to a neat way of solving linear Diophantine equations. 
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Example 1.4 


Find the solution of the linear Diophantine equation 13x + 59y = 10 in 
which x takes its smallest positive value. 


We begin by finding the FCF of 5, 
59=4x13+7 
13=1x7+6 
T=1x6+1 
6=6x1 
59 45 5+4 9 
— = 1 i —, —, OCU = — 
So 3 [4,1,1,6], and its convergents are os rs ae and 
6x9+5 59 


a> = 79: From the final two convergents, (as Theorem 1.3 
Gx. T3 
guaranteed would occur), we observe that 
59x 2—9x13=1. 
Multiplying through by 10 and rearranging gives 
13 x (—90) + 59 x 20 = 10, 


showing that x = —90, y = 20 is one solution of 132 + 59y = 10. The general 
solution is therefore 


r=-904+59t, y=20-13, teZ. 


For the smallest positive value that x can have in a solution take t = 2, to 
give the solution z = 28, y = —6. $ 


Problem 1.8 


48 
By determining the convergents of T9’ find the general solution of the linear 


Diophantine equation 48x + 19y = 4. 


Returning to Theorem 1.3, the alternative formulations of properties (a) and 
(b) give important information about the convergents of an FCF. 
Property (a) tells about the difference between successive convergents: 


(=i 
Cy = C= = Leia eat E E A (*) 
dklk-1 


Now the denominators qk satisfy 
Mm S92< 93 <-++<Q. 


Therefore the right-hand side of equation (x) becomes smaller as k increases. 
That is, the convergents are coming closer together. 


Property (b), on the other hand, tells about the difference between 
convergents two apart: 
zta 
Beh og = FORO — pA a tae (+) 
Gk k—2 
We shall use both these facts in the proof of the following theorem 
concerning the relative sizes of convergents. For ease of notation we refer to 
the convergents C2, C4, C6, ..., which have even subscripts, as the even 
convergents whilst C1, C3, C5, ... are called the odd convergents. 


We could equally well have used 


the convergents of z but taking 


the larger number on the top is 
more efficient. 


Remember that Ck = Pee, 
dk 


We ask you to verify this in 
Problem 1.10 
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Theorem 1.4 The relative sizes of convergents 


For any FCF the odd convergents form a strictly increasing sequence 


C1 < C3 < Cs <---, while the even convergents form a strictly 
decreasing sequence C2 > C4 > Ce >~.. 


Every even convergent is greater than every odd convergent. 


We can illustrate pictorially what the theorem is claiming by imagining a 
plot of the points C, against r. 


ee a a x Cn 


Cı 
Figure 1.3 The convergents of an FCF 


The final convergent Cn, which is the value of the FCF, may be either an 
even or an odd convergent. The last part of the theorem tells us that 
whichever it is, Cn is at least as large as every odd convergent and no larger 
than any even convergent. 


Proof of Theorem 1.4 
As ax > 0 (for k > 2) and qp > 0 (for k > 1) we note that, for k > 3, 
(Dla; 
dkdk-2 


has the same sign as (—1)*~! and so will be positive when k is odd and 
negative when k is even. 


Putting k = 2r + 1 in equation (**) produces 

Catt — Cor-1 >0, forr>1, 
which confirms that the odd convergents are strictly increasing. 
Putting k = 2r in the same equation gives 

Cor — Cop_2 < 0, forr>1, 


which confirms that the even convergents are strictly decreasing. 
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It remains to show that every odd convergent is less than every even 
convergent, and for that we use equation (*) which tells us that Ck — Ck—1 
takes the same sign as (—1)*, for k = 2, 3, ..., n. Putting k = 2r gives 


Cor — Co 0, for r >, 
and putting k = 2r + 1 gives 
Cor41 — Corp <0, forr>1. 


Hence the even convergent C2, is larger than each of the adjacent odd 
convergents Co,_; and C2,41. 


Now consider any even convergent C2, and any odd convergent C241. 


If s > t then Co, > Cos41 > Cor41, the second inequality utilizing the fact 
that the odd convergents increase. 


If s < t then Co, > Ca > C2441, the first inequality utilizing the fact that 
the even convergents decrease. 


In each case C2, > C2441, as required. 


Problem 1.9 


Determine the convergents of [2, 1, 1, 2, 1, 1, 2] and confirm that they 
satisfy Theorem 1.4. 


Problem 1.10 


Show that, for any FCF, the denominators q, are increasing, and are strictly 
increasing with the possible exception of q1 = q2. Show further that, for any 


FCF, gq, > k for all k > 5. 
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2 INFINITE CONTINUED FRACTIONS 


2.1 The convergence of a continued fraction 
In the introduction to the unit we built up a sequence of FCFs 
ah AEE 
whose rational values give better and better approximations to the golden 


+v5 


ratio is irrational, and so no matter how far we go along this sequence of 
FCFs we can never reach the golden ratio itself. What is of interest is the 
limit of this sequence of rationals. This suggests that we might explore the 
idea of the unending continued fraction [1,1,1,...]. 


ratio, 


. However, the value of any FCF is rational, but the golden 


This is our first example of a simple infinite continued fraction. 


Definition 2.1 Infinite Continued Fractions. 


Let a; be any real number and let a2, a3, ..., an, ..., be any positive 
real numbers. The expression 


1 


a, + 
a2 + 
a3 + 


re 
Ol = 


is an infinite continued fraction. We write this continued fraction in 
terms of its partial quotients as [a1, a2, a3,...]. Its convergents are the 
finite continued fractions 


C1 = [ai], C2 = [a1, a2], C3 = [a1, a2, a3],.... 


When all the partial quotients are integers, the continued fraction is 
said to be simple. We shall use the abbreviation ICF for a simple 
infinite continued fraction. 


In this subsection we shall explore basic properties of infinite continued 
fractions. Our ultimate interest will lie with simple infinite continued 
fractions although, as in the previous section, we shall have occasion to allow 
the partial quotients to be non-integers for the purpose of providing general 
proofs. 


What exactly does an infinite continued fraction represent? Whereas 
calculating the value of a finite continued fraction involves only the normal 
operations of arithmetic, evaluation of an ICF is an infinite process and as 
such requires a different approach. There is a strong resemblance to the 
convergence question for infinite series. The sum of the infinite series 


co 
Ņ ap = t1 +22 +23 +- 


r=1 
is defined to be the limit of the sequence 
Cy, Cito, C1 Fee es, 22.5 rtr teers ain, 


of partial sums, provided this limit exists. For ICFs the convergents fulfil 
the role of the partial sums and so, by analogy, we define the value of an 
ICF to be the limit (if it exists) of the sequence {Cn} of its convergents. 
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We use the notation {Cn} to 
represent the infinite sequence 


(Gul OF Cpo: 


Definition 2.2 Value of an Infinite Continued Fraction 


The value of the infinite continued fraction [a1, a2, a3,...] is defined to 
be the limit of its sequence of convergents {Cp}. 


Notice that we did not add the condition ‘provided this limit exists’. This is 
because there is one big difference between infinite series and ICFs: whereas 
some infinite series converge and some do not, it transpires that for every 
ICF, its sequence of convergents has a limit. 


Theorem 2.1 Sequences of convergents converge 


For any ICF [a1, a2, a3, ...] its sequence of convergents {Cn} tends to a 
limit. 


Proof of Theorem 2.1 
For any n > 1, the the first n convergents of the ICF [a1, a2, a3,...] are 
Ci, C2, ..., Cn—1, Cn. These are precisely the convergents of the FCF 


Cn = [a1, a2, a3,..., an]. What is more, Cn = Bs. where pn and qn are the 
numerators and denominators as determined by the formulae in 
Theorem 1.2. 
Now Theorem 1.4 tells us that 
Cy < C3 <--+ < Congr <--> < On < +++ < Cop <--- A < Cd, 


where 2k + 1 <n. As such a chain of inequalities holds for any value of n, 
no matter how large, we see that the odd convergents of the ICF form a 
sequence C1, C3, Cs, ... which is strictly increasing and bounded above 
by C2. A well-known theorem of analysis tells us that a strictly increasing 
sequence of real numbers which is bounded above is convergent, so 

Ci, C3, C5, ... has a limit, say a. 


In the same way the even convergents form a strictly decreasing sequence 


C2, C4, Ce, ... which is bounded below by C1, and so it has a limit, say b. 
iG 
e 
N 
w 
* 
N 
s E j Ch 
Me steers ofy Cs b 
_--+---- =. =a 
Pn ess ae Cs 
- 2 E a 
Pie C3 
5 a 
Z 
d 
ra 
£. 
r 
y: 
LA 
€ 
vA 

é 

Cı 

Figure 2.1 


We shall use the term ‘converges’ 
as an alternative to ‘tends to a 
limit’. So the following theorem 
tells us that the sequence of 
convergents of every ICF converges! 


This result usually goes by the 
name of The Monotone 
Convergence Theorem. Here, we 
are going to assume this result is 
known, along with some other 
results on the convergence of 
sequences of real numbers. 


21 


For the whole sequence C1, C2, C3, C4, ... to converge it remains to show 
that a = b. As every even convergent is greater than every odd convergent, a 
and b both lie between C2,_1 and Con, for each positive integer n. Now 


1 
Con — Con-1 = ————._ by property (a) of Theorem 1.3, 
G2nQ2n—-1 


1 
<< (2n)(2n — 1)’ since dk = k for large k Problem 1.10. 
Therefore 
1 
(2n)(2n — 1)’ 
and, as the term on the right tends to 0 as n increases without bound, we 
conclude that a = b. a 


la ae b| < Con i Con-1 = 


The above proof does more than confirm that the limit exists; it contains an 
estimate for the accuracy of the convergent Cn as an approximation to the 
actual limit. We know that the limit lies somewhere between any pair of 
consecutive convergents. If we couple this with the information given by 
Theorem 1.3 property (a), which gives a bound for how far apart successive 
convergents can be, we have the following immediate consequence. 


Corollary to Theorem 2.1 When we talk of the convergents 

of x we mean the convergents of 

If {Cn} is the sequence of convergents of x then, for each n > 1: the ICF of zx. 

(a) x lies between Cn and Chn+1; 
1 1 


5 ia = 
QnQn+1 


dn 


Punti _ Dn 
dn+1 dn 


T 


Qn Qn4+1 j 


Problem 2.1 


Estimate the maximum possible error in taking the convergent Cs as an 
approximation to the value of [1,1,1,1,...]. What is the maximum error in 
Ce? The actual value of this ICF is the golden ratio. Taking the latter to be 
1.61803 to five decimal places, what are the actual errors in Cs and Cg? 


2.2 Periodic continued fractions 


We now know that any ICF has a value, but how do we go about finding that 
value? One thing we can do is determine more and more of its convergents 
in the knowledge that we are obtaining better and better rational 
approximations to its value. More about that later, but for the time being 
there are some ICFs whose value we can find without too much difficulty. 


Example 2.1 


Find the value of the ICF [1,3,2,1,3,2,1,3,2,...], where the sequence 1,3, 2 
of partial quotients is repeated indefinitely. 


Knowing that the ICF converges, let its value be x. Then notice that 


1+— ——_ = 1+ —__,_ = 1+ 
D 3+ 
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and, as this continued fraction has value z, 

g= iia r]: 
The convergents of [1, 3,2, x] are, using the formulae for the numerators and 
denominators, 

E A <9 9x +4 

TRTA 
and, as the final convergent is x itself, we have the equation 

97 +4 

ces. TE 

Clearly x > a, = 1 and so we want the positive root of this equation. Using 


a ‘ 


that is 7x? — 6x — 4 = 0. 


the formula for the roots of a quadratic equation we obtain z = 


Problem 2.2 


Consider the ICF [3, 1,3,3, 1,3,3,1,3,...], where the sequence of partial 
quotients 3, 1, 3, recurs indefinitely. Determine its value x by writing 


1 


1 1 
trey 


3+- 
T 


r=3+ 


and solving the resulting quadratic equation. 


Problem 2.3 


Find the value of the ICF [2,2,2,2,2,...]. From this result write down the 


i 
ICF for v2 and hence find a rational number z which lies within —— 


1000 
of v2. 


In Example 2.1, where we effectively wrote 
t=1+ 

Pa 

£ 


we were taking some liberties with properties of limits. To give the required 
argument correctly we should really be dealing with finite expressions, and 
subsequently taking limits in the knowledge that they do exist. That is, 
instead of working with the ICF whose value is x we should deal with its nth 
convergent, Cn. From our work on FCFs we can certainly write 


Crh=1+ 5 , forn>4, 
3+ z T 
F fs 
which leads (as above) to 
9Cn-3 + 4 
ek ES 


Now as n increases without bound we know that Cn converges to x, and 
consequently C,3 converges to x. Taking limits of both sides gives 


9r +4 
Ss s 
Tz +3 
which is the situation we reached in Example 2.1. 
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The property which we have been using is the following. 


Infinite continued fraction identity 


1 


lai, a2, 43, .. | Sis eS ee 
[a2, a3, a4,...] 


This follows since we have the corresponding result for FCFs 
aT 2 See 
laz, AZ3,+++5 an] 


and, as n increases without bound, the FCFs on each side of this equation 
converge to the required ICFs. 


EA =a, + 


So far, the only ICFs which we have been able to evaluate are ones in which 
a block of partial quotients repeat indefinitely. We have a name and 
notation for such ICFs. 


Definition 2.3 Periodic Continued Fractions 


An ICF which consists of a block of partial quotients which repeat 
indefinitely is said to be purely periodic, and we write it as 
[(a1, a2, a3, . - -,an)]| enclosing the repeating block in angular brackets. 


We refer to the repeating block as the cycle of the ICF. 


An ICF of the form [b1, b2, .. . , bm, (a1, a2,- - - , an}], in which a 
repeating block occurs after a finite number of partial quotients, is said 
to be periodic. 


For example, 
BASAS A ISKA 
is a purely periodic ICF, and 
Monne eee o ae | ER] 


is periodic. Periodic continued fractions may be written in different forms 
using the angle bracket notation. For example, 


[(3, 4)] = [3, (4, 3) = [3, 4, (3, 4)]. 


The evaluations that we have carried out illustrate a general strategy for 
finding the value of a purely periodic continued fraction. To summarize, 
consider |(a1, a2, a3, ...,an)}| with value x. The key step is the observation 
that x can be written as the finite (but not simple) continued fraction 


T= le 05, 43,0550 52) 


Adopting the usual notation for the convergents, the final two convergents of h 
Pn-1 


and Pas and can be determined by the methods 
Qn-1 dn 

of the previous section. The final convergent of [a1, a2, a3, . . - , an, £] is then 

TPn T Pn-1 


Ldn + In—-1 


[a1, a2, a3, . . . , an] are 


, and equating this to x yields the quadratic equation 


gar? $ (Qn—1 = Dn)& — Pn—-1 = 0. 


As the first partial quotient a; reappears infinitely often as a partial 
quotient it has to satisfy a; > 1, and so the root we seek is certainly 


positive. In fact this quadratic has one positive and one negative root since 
Pn-1 


the product of its roots is equal to — which is negative since pn_; and If az? + br +c = 0 has roots a and 
dn 


c b 
qn are each positive. B then aß = = anda+6= ee 
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It is only a short step from the evaluation of purely periodic continued 
fractions to the evaluation of any periodic continued fraction, as illustrated 
by the following example. 


Example 2.2 

Evaluate the periodic continued fraction [—2, 1, (3, 4)]. 

Let y = [-2, 1, (3, 4)] and x = [(3, 4)], so that y = [—2, 1, z]. 

We evaluate the purely periodic x by noting that £ = [3, 4, x], the latter ICF 


3-1 1 3 
having convergents Da and a Equating the final convergent to x 
gives the quadratic 
4a? — 12r —3 =0, 
2v3 
and this has positive root z = = 
—2 — -x-2 
Turn now to y = [—2,1,2]. It has convergents T and 1” and so 
3+4+2V3 
<< ae ee eee 
eet S98 - 8423 
—— EA 
-27223 5-273 
5+2V3 5-2/3 
_ —23+4V3 + 
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The method illustrated in Example 2.2 can be used to evaluate any periodic 
ICF. 


Problem 2.4 
Evaluate the following ICFs. 


(a) [2,(1,1,1,4)] (b) [2,2,(3)] 


Note carefully the final step in the 
calculation in this example. In 
a +byn 
dyn 
may be simplified by multiplying 
top and bottom by c — dyn and 
appealing to the identity 
(c + dyn)(c — dyn) = c? — ng. 
The net result is that all 


square-roots will be removed from 
the resulting denominator. 


general the expression 
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3 THE CONTINUED FRACTION 
ALGORITHM 


3.1 The continued fraction of an irrational 


Having seen that every finite simple continued fraction has a rational value 
and, conversely, that every rational number can be represented by a finite 
simple continued fraction, the following result should come as no surprise. 


Theorem 3.1 Irrationality of ICFs 


The value of any infinite simple continued fraction is an irrational 
number. 


You may have thought that this theorem must follow immediately, but we 
must exclude the possibility that a rational number could be represented by 
both a finite and an infinite simple continued fraction. 


Proof of Theorem 3.1 


Let x = [a1, a2, a3, . . -| be an ICF. For each n > 1, x lies between the 
consecutive convergents Cn = Pn and Cati g Nei and so 
Qn Qn+1 
Pn 1 
r— — : 
dn Qn Qn4+1 


7 

Now suppose that x is rational; say x = —, where r and s are integers with 
8 

s>0. Then 


r TQn — sS 1 
o< |" — Pa| — Irae Pnl < : 
s dn Sdn QnQn+1 
so that 
0 < |ran — SPn| < 
Qn+1 


But the denominators qn form a strictly increasing sequence of integers and 
so we can choose n so that qn+1 > s. But then |rqn — spn| would be an 
integer lying strictly between 0 and 1, which is impossible. The conclusion 
from this is that no such integers r and s can exist; in other words, x must 
be irrational. a 


We turn now to the converse problem. Given an irrational number can it be 
expressed as a simple infinite continued fraction? In answering this question 
in the affirmative we shall discover more: the simple continued fraction of an 
irrational number is unique. Our approach will be to demonstrate a method, 
called the Continued Fraction Algorithm, for determining an ICF for any 
given irrational number. Of course we cannot determine all its partial 
quotients unless there is some recurring pattern such as occurs with periodic 
ICFs. What the algorithm does is enable us to determine as many partial 
quotients as we please. The situation is analogous to the problem of finding 
the decimal expansion of 7; the full expansion of m cannot be given, but 7 
can be determined to any required number of decimal places. 


Given an irrational number zı we seek integers a1, a2, ... so that 
r= [a1, a2, os ale 


Knowing that zı lies strictly between any successive pair of its convergents, 
1 
it has to lie between C1 = a; and C2 = a, Be —. In particular, as a2 > 1, 7 


lies strictly between the integers a; and a, + T It follows that a; = int(x1), 


26 


See the Corollary to Theorem 2.1. 


Notice that 

1 = Pn 

S dn 
for otherwise the continued fraction 
would be finite. 


Remember int(x) denotes the 
integer part of x. 


and we note, in particular, that a is determined uniquely from 2 alone. 
Now let us write 


zı = int(x1) + frac(z,), 
where frac(x,) denotes the fractional part of xı. Then 


1 > 1 5 1 
Pea 7 E int(x1) + a int (a) 


7 a2 + 
Be — Ce Na —— 


+ 2 ee 
[a2, a3, @4,...| 


and so 

1 1 í oda 
SSS ESS The fractional part of an irrational 
zı —int(xı) frac(zı) number lies strictly between 0 and 
1, and so its reciprocal is greater 
then z2 = [a2, a3, a4, . . .| and x2 > 1, as it is an than 1. 


[a2, a3, a4, Sy | = 


If we write ry = : 
— frac(x1) 


irrational number greater than the positive integer az. Repetition of the 
above argument now leads to ag = int(x2), and so a2 is determined uniquely 
from x2 and therefore from zı. Continuing in this way, writing 23 = ———— 
frac(x2) 
we have z3 = |as, a4, a5, ...] and a3 = int(r3), etc. This is the Continued 
Fraction Algorithm which is summed up in the flowchart of Figure 3.1. 


Irrational xı 


Express this number as 


integer part + fractional part 


Integer part is the next Calculate the reciprocal 
partial quotient = of the fractional part 


Figure 3.1 The Continued Fraction Algorithm 


Ze 


3.2 Using the Continued Fraction Algorithm 


Before pursuing the theory let us see the algorithm in action. 


Example 3.1 


Determine the continued fraction of V6. 


Suppose that zı = V6 = [a1, a2, a3,...] and we are to use the algorithm to 
determine these partial quotients. Since 2? < 6 < 3?, V6 lies between 
2 and 3. Therefore int(/6) = 2 and V6 = 2 + (V6 — 2) expresses it as 

1 


integer part, 2, plus fractional part (v6 — 2). So a; = 2 and x2 = v r 


where £2 = [a2,a3,...]. 
To find az we require the integer part of x2. 


1 1 V6+2_ V6+2 


T2 = —=— = x = F Note the trick of multiplying top 
v6-2 V6-2 V6+2 2 and bottom by V6 + 2 to remove 
Now, as V6 lies between 2 and 3, v6 + 2 lies between 4 and 5, and so the square roots from the denominator. 
2 
integer part of is 2. Therefore 
6-2 
anpi (*) 


expresses X2 as its integer part plus its fractional part. So ag = 2 and 
2 


T3 


= e E 
We continue by finding the integer and fractional parts of x3: 
2 2 V6+2 
To = x = V64+2=44 v6 =~ 2). ** The integer part of V6 + 2 is 4. 
‘= eee eae ( ) (**) ger p 
1 
Solas —4.and 4. — oe But now 24 = £2 and we are in a cycle. We 


repeat calculation (*), establishing a4 = 2, then repeat (**), establishing 
as = 4, then back to (*) and so on. The sequence of partial quotients 
repeats the cycle 2, 4 indefinitely from here onwards and we conclude that 


v6 = (2, (2, 4)]. + 


Our second example is similar, but this time we shall give just the 
calculations involved without the commentary. 


Example 3.2 
Determine the continued fraction of : ae 
7+ 2 V2-1 2 
= =4 i. =4 = 
ry 2 ae 2 , ay jee Vac] 
2 V2+1 1 
T2 = x = 2f/24+2=44 (2V2-2); ag=4, 23 = —~— 
t= l  2v242 _ vV2+1 yl, a DE z 
3 2/2 —2 W2+2 2 2 ’ 3 ’ 4 vii 
L4 = LQ, and sO a4 = G2, G5 = a3, Ag = a4, etc. 
Therefore 
7+ V2 
v= 14, (4,1) + 
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Problem 3.1 


Use the Continued Fraction Algorithm to determine ICFs for the following 
irrational numbers. 


(a) V3 (b) 24+ V7 


Having seen illustrations of the Continued Fraction Algorithm in use, we 
must now do a little tidying up. We have seen that the algorithm will always 
produce an ICF, but can we be certain that the value of this ICF is equal to 
the irrational number from which we started? There is also a question of 
uniqueness to be resolved; is the ICF of an irrational number unique? 


Theorem 3.2 Uniqueness of the ICF 


Given any irrational number x, the ICF determined from x by the 
Continued Fraction Algorithm converges to x, and no other ICF 
converges to T. 


Proof of Theorem 3.2 


From the irrational number x the Continued Fraction Algorithm produces 
a sequence of (non-simple) finite continued fractions as follows: 


g= 21 = (a1, L2| = [a1, a2, £3] = [a1, a2, 03, 24] = --- = [A1, G2, a3,.--, An, Ln41], 


where, for all n > 1, an = int(x,,) and the step from one stage to the next is 
accomplished by 


Ly — Gy + 7 
Tn+1 


To show that the ICF which is being built up in this way has limit x we 
must show that the sequence of finite continued fractions 


[a1], [a1, ag], paths (2102; 2,01; 
converges to T. 


Using the standard notation for the numerators and denominators, 


[ain Go,. 20; ax| — a we have for n > 2 
k 
Pn l ul 
r——|< roy By the Corollary to Theorem 2.1 
dn QnQn+1 


As the denominator sequence {qn} is, for n > 2, a strictly increasing 


tends to 0, confirming that =~ can be 
QnQn+1 dn 
made as close as we please to x by choosing n sufficiently large. Hence the 


ICF determined from x by the algorithm does indeed converge to x. 


sequence of positive integers, 


It remains to show that the ICF of an irrational number is unique. To that 
end suppose that the ICFs [a1, a2, a3, ...] and [b1, be, bs,...] have the same 
value. We shall use the Second Principle of Mathematical Induction to show 
that an = bn for all n > 1, and so the two ICFs are identical. 


First we observe that if [a1, a2, a3,...| = [b1, b2, b3,...] then a; = bı since 
both are equal to the integer part of the common value. This gives the basis 
for the induction. 


For the induction step suppose that, for some k > 1, 
ay Eos a2 =b, cosy Ak = pp; 
Our goal is to deduce that a441- = bkr- Now, 


[a1, a2, a3, . x | at [a1, Q2,. ++ ,Qk, [ak+1, h+42, =e |] 
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and similarly 
[b1, ba, bs, ---] = [b1, b2,- 


As these two have the same value, and have the same first k partial 
quotients, it follows that 


«4 Oks [Ont 5 O25 af: 


[Qx+1,0%42,---] = [be41,0n+42,.-.]. 


But now ax+41 = bk+1 as each is equal to the integer part of the value of this 
ICF, and the induction step is complete. Hence the two ICFs having the 
same value are the same. a 


3.3 The continued fraction of yn 


Let us take stock of where we have reached with simple continued fractions. 

For FCFs, the value of any FCF is rational and conversely every rational can 
be represented as an FCF. However, the representation is not unique; indeed, 
each rational number can be represented as an FCF in exactly two ways. 


Turning to ICFs, the value of any ICF is irrational, and conversely every 
irrational number can be represented as a ICF. But this time the 
representation is unique; every irrational number has a unique ICF. 


The only ICFs which we have been able to evaluate so far have been 
periodic ones, and these have always turned out to involve square roots in 
some way. We are not going to digress to prove it here, but it turns out that 
the values of periodic continued fractions are precisely the numbers of the 
form r + syn, where r and s are rational and n is a non-square positive 
integer. That is, each periodic ICF has value equal to one of these numbers 
and, conversely, each such number has a periodic ICF. For the present, let us 
investigate further the ICF of certain square roots. 


Example 3.3 
For each positive integer n determine the ICF for Vn? + 1. 


The integer part of Vn? + 1 is n, and so we write 
Vn?+1=n+ (vn?+1-n). 


This gives aj = n. We look next at the reciprocal of Vn? + 1 — n, which is 
simplified by multiplying top and bottom by Vn? +1-+n. 


n =m 


So ag = 2n, and the same fractional part vn? + 1 — n has arisen. This 
means that we are immediately in a cycle and so vn? + 1 = [n, (2n)]. 


For example, for n = 1, 2, 3 and 4 we obtain, respectively, the following. 
V2 = [1, (2) 
v5 = [2, (4)] 
v10 = [3, (6)] 
V17 = [4, (8)] + 
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This proof is formalizing the 
repeated use of the infinite 
continued fraction identity. From 
[a1, a2, a3, . =] = [b1, b2, b3, oe =| we 
deduce that a; = bı and 

[a2, a3, a4, . jl = [b2, b3, b4, oe als 
But then az = be and 

[a3, a4, -$ J = [bs, Das .] etc. 


These numbers are called the 
quadratic irrationals and arise as 
roots of quadratic equations 

az? + bz +c = 0, where a, b and c 
are integers. 


Problem 3.2 
Determine the continued fraction of Vn? + 2 for any integer n > 1. 


Hence find ICFs for, /3, V6, V11 and v18. 


By now you should be feeling confident that you could use the Continued 
Fraction Algorithm to determine the ICF of yn, for any positive integer n. 
In the table below we collate some of the ICFs of y/n which we have seen 
and give several more. 


Table 3.1 Simple continued fractions for yn 


v2 = [1, (2)] 713 = [3, (1, 1,1,1,6)] V23 = (4, (1,3, 1,8)] 
v3 = (1, (1, 2)] v14 = [3, (1,2, 1,6)] v24 = [4, (1,8)] 
v5 = [2, (4)] v15 = [3, (1,6)] v26 = [5, (10)] 
v6 = (2, (2,4)] v17 = [4, (8)] v27 = [5, (5, 10)] 
VT = [2, (1,1, 1,4)] v18 = [4, (4, 8)] V28 = [5, (3, 2,3, 10)] 
V8 = (2, (1,4)] V19 = [4, (2,1,3,1,2,8)] v29 = [5, (2, 1,1, 2, 10)] 
V11 = [8, (3, 6)] Ji = (4, (1,1,2,1,1,8)] V31 = [5, (1,1,3,5, 
v12 = [3, (2, 6)] v22 = [4, (1, 2,4, 2,1, 8)] v32 = [5, (1,1,1, J 
V46 = (6, (13,112.62 1 1 3,1, 12)] 
V61 = (7, (1,4,3, 1,2, 2, 1,3, 4, 1, 14)] 
vos = 19, (1,2,9.1,3-5,1,8, 4,5, 9, 1,0, 2,509) 
497 = (9, (1,5, 1,1,1,1,1,1,5, 1, 18)] 
VIEI = (12, (3,2, 7, 1,3,4, 1, 1, 1, 11, 1, 1, 1,4,3, 1, 7, 2,3, 24)] 

166 = [12, (1,7, 1,1,1,2,4, 1,3, 2, 12, 2,3,1,4,2,1,1,1,7, 1, 24)] 

Problem 3.3 


Examine Table 3.1. 

(a) What common features do you observe in the continued fractions? 

(b) Can you draw any conclusions concerning the continued fraction of 
Jn + int(/n)? 

(c) Conjecture the general form of the continued fraction of y/n, for n a 
positive non-square integer. 


From our solution to Problem 3.3 we are led to believe that, for any 
non-square positive integer n, the ICF of yn takes the form 
Vn = [a1, (a2, a3,...,@3, a2, 2a1)], 
where a; is the integer part of y/n. Here the periodic part may consist of: 
e just (2a;) as in, for example, V2 = [1, (2)] and V17 = [4, (8)]; 
e an odd number of partial quotients as in V13 = [3, (1,1, 1,1, 6)]; 
e an even number of partial quotients as in V7 = [2, (1, 1,1, 4)]. 
It transpires that the general form is, indeed, as claimed. However proof of 


this would take us deep into the study of more properties of continued 
fractions, which we do not have space for in this course. 


æ oO 
w 
ran 
en 
= 
© 
e 
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4 RATIONAL APPROXIMATIONS 


4.1 Approximating irrationals by rationals 


We have seen that an irrational number has a unique ICF and that the 
convergents provide a sequence of better and better rational approximations 
to the irrational number. We even have a way of estimating how good these 
rational approximations are. In this section we shall go further and show 
that continued fractions provide what are, in a certain sense, best rational 
approximations. 


Given an irrational number z, we can always improve on any rational 
À A a š 4 2 
approximation D to x. In fact there are infinitely many rational numbers 


lying between ; and x. Consequently there is no such thing as the best 


rational approximation to x. However, if we impose restrictions on the size 
of the denominator we can ask about best approximations. For example, of 
all rational numbers whose (positive) denominator does not exceed 50, the 


2 
best rational approximation to 7 is 7 We have to allow the denominator 


179 
to climb to 57 before we find the better approximation ET to T. 


What we shall show in this section is that if Pr is a convergent of x then 
In 


a 
there is no rational number ; with 0 < b < qn which approximates x better 


than ee: and in this sense the convergents of x are best rational 
q 


approximations to it. But before becoming involved in the analysis we give 
an elegant geometric interpretation of the convergents to a continued fraction 
of an irrational number which was first given by Klein in 1895. It illustrates 
precisely the sense in which convergents are best rational approximations. 


Let a > 0 be irrational and consider the part of the line y = az in the first 
quadrant of the Cartesian plane as illustrated in Figure 4.1. This line cannot 
pass through any of the lattice points (i.e points whose coordinates are both 
integers) except (0,0), for if it did œ would be rational. Imagine that pins 
are set upright at each of the lattice points and a fine string is placed along 
the line with one end fixed to the line at an infinitely remote point, and the 
other end at the origin. Keeping the string tight we move the end at the 
origin towards the right. It will catch on a series of pins below the line. In 
fact, the pins on which it catches are those at lattice points with coordinates 
(q1,P1), (43, p3), (q5, Ps), --- corresponding to the odd convergents 

pı Ps Ps 

d 93° 45° 
of a. Similarly, if we move the string from the origin to the left it catches on 
pegs at (q2, p2), (q4, pa), (q6, Pe), -.., corresponding to the even convergents 

p2 pa Pe 

qo’ qa’ q 
of a. 
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Remember that when we write i 


for a rational number we assume 
that a, b are integers with b 
positive and gcd(a, b) = 1. 


If y = ax passes through (m,n) 
then a = £. 
m 


Recall that 
Ci < C3 < C5. Ae Cases? O < G, 


and notice that C; = p: give the 


gradients of the lines joining the 
origin (0,0) to the points (qi, pi). 


5 


Figure 4.1 Convergents as best approximations 
to the golden ratio 


1 5 
Figure 4.1 illustrates all this for the case a = EALA the golden ratio, 
I erro Dg 1S 
h in =, =, =, =, 2, —,.-. . Thed ine i 
whose convergents begin Ea 8” e dashed line is the 


portion of y = az in the first quadrant. Moving to the right it catches on 
pins with coordinates (1, 1), (2, 3), (5,8),..., whilst moving left it catches on 
the pins at (1,2), (3,5), (8,13), and so on. 


Let us return to the question of best rational approximations. First we must 
prove, formally, a result which we have seen illustrated in the many 
sequences of convergents that we have calculated; each convergent is a better 
approximation to the irrational number than is its predecessor. The theorem 
which follows includes this result, but also gives some bounds for the 
accuracy of the convergents. 


Theorem 4.1 Accuracy of convergents 


If se = s ... is the sequence of convergents of x then, for all k > 1, 
di. g7 493 


< 


k- ood eget Sears Le ae 
dk+1|  Qk+19k+2 — 20k Qk+1 
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Proof of Theorem 4.1 


Let x have ICF [a;, a2, a3,...]. The Continued Fraction Algorithm yields the 
system of equations 3 


x = [z2] = [a1, a223] = --- = [a102i an Baral =-°- 
and 
oo Pa = [@1,@2,... An Enri = Pn 
n dn 
= Tn+1Pn + Pn-1 = Pn f by Theorem 12 
Tn+419n + Qn-1 Qn 
= Pn—-19n — Pndn-1 
Qn(Zn+19n F Gn—1) 
1 
= ——, by Theorem 1.3 property(a). 
Qn(n4+19n + Qn—1)’ 4 ( ) 
Now Tayi = [€n+41,4n42,4n43,---] and 80 Qn41 < Enpi < ang +1. 
Therefore 
Pn 1 1 
T — ’ 
Qn | — Qn(@n419n+4n—-1) Gn Qn41 
which gives the left-hand inequality when n = k + 1. But we also have 
1 
Qn dn((an+1 + 1)dn + Gn—1) 
a 1 n 1 
Qn(On414n ae Qn—1) + g Qn(Qn+1 oe In) 
i 1 
> — |, = Remember that {qn}, for n > 2, is 
Qn(Qn+1 + 4n+1) — 24nQn+1 a strictly increasing sequence. 


which gives the right-hand inequality by taking n = k. 
For the middle inequality note that 


Gk+2 = Ak+29k41 + qk = qk + Ok = 24k 


and so 
1 1 


dk+14k+2 — 2qk+1qk 


The left-hand inequality in Theorem 4.1 gives an indication of how close 
each convergent is to the true value whilst the right-hand inequality gives a 
bound which limits the accuracy. We can collate the two facts for each 
convergent as follows. 


Corollary to Theorem 4.1 


For all k > 1 the convergent Pk of x satisfies 
dk 


1 1 
<|e- 2< 


2dkqk+1 dk| © dkqk+1 


Problem 4.1 


The ICF of V5 is [2, (4)]. Estimate the accuracy of the third convergent as 
an approximation to v5. 
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4.2 Best approximations 


The discovery that the ICFs of y/n can be described exactly in terms of 
finite recurring cycles is an interesting fact, because the decimal 
representations of these numbers (when n is non-square) have no discernible 
pattern. Of course ICFs may have other patterns besides recurring blocks 
which permit them to be described exactly. Amazingly, one of the 
best-known of all irrational numbers, e, the base of the natural logarithms, 
has such a pattern to its ICF. Euler discovered in 1737 that 


e = [2,\1,2,2,1,4)1)1;6,1,4,8,1,1,1043.....1, 


where the even positive integers 2, 4, 6, ... occur in order as every third 
partial quotient, the other partial quotients all being 1 with the exception of 
ay = Ds 

Example 4.1 


Determine the first convergent of e which is accurate to within 0.0001 of the 
true value of e. 


The convergents of e are 


2 3 8 11 19 87 106 193 1264 
1 Loe ae ee T; a 


11 
Using the result of the Corollary we know that the convergent E lies 


1 19 1 
within Tae of e, 7 lies within PET of set so on. The first convergent 
guaranteed to be accurate to 0.0001 of e is 7 which lies within 
1 


106 
= 0.0000303... of e. The previous convergent 30 does not have The values to six decimal places 


71 x 465 
the required accuracy because, using the left-hand inequality in the R eee 
e=z. 
Corollary, the error in 39 8 at least > = 2.718310 
106 
iPibienit42: $$ ee ee ee 
The ICF of V2 is [1, (2)]. Determine the first convergent which is accurate 
to within 0.0001 of the true value of V2. 
. 193.. AA 

We now appreciate that the convergent —— is a very good approximation 
to e, but as yet we have made little progress towards establishing the sense 
in which it is ‘best’. We shall put this right immediately by showing that of 

r ? : 7 193 
all rational numbers with denominator not exceeding 71, the convergent 71 
is the best approximation. The general result from which we can draw this 
as a conclusion follows quickly from the slightly stronger result of 
Theorem 4.2 below. 
For completeness we have given proofs of all the remaining results in this 
section but you should regard these proofs as optional reading. 

Theorem 4.2 
Let 2" be the nth convergent of the irrational number z, and let = be Remember = is in lowest terms. 


dn 
any rational number with 0 < b < qn+1. Then, for n > 1, 


lanz — Pn| < |ba — al, 


the equality holding only if a = pn and b = qn. 
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Proof of Theorem 4.2 


a 
Let 3 be a rational number (in its lowest terms) with b < qn+1, and suppose 
it is not the case that a = pn and b = qn, when the required equality 
certainly holds. 

Consider the system of equations (with r, s real for the moment): 


Q = TPn + SPn+1; 
b = Tn + 8Qn41- 


Multiplying the first equation by qn and the second by pn and subtracting 
leads to 


aqn — DPn = 8(Pn+19n = PnQn+1) 


which, after applying Theorem 1.3 property (a) to the bracketed term on the 
right-hand side, gives 


s = (—1)”+! (aqn — bpn). 


Similarly 


Mo (—1)"**(6pn4i — 0Qn+41)- 


Note that r and s are therefore integers, and neither of them is 0 because 
the assumption that r = 0 leads to agn41 = bpn41, whereupon Euclid’s 
Lemma forces qn+1 to divide b, since gcd(Pn+1, qn+1) = 1, and this 


contradicts 0 < b < qn+1. Also if s = 0 we have i = D which is the case of 
equality which has been excluded. e 


Now since 0 < b = rqn + Sqn+1 < qn+1, the non-zero integers r and s must 


have opposite sign. Moreover x — Pn etl e 
dn+1 


because the convergents are eee greater than and then less than 
the irrational limit x. It follows that r(qn£ — pn) and s(qn+1£ — Pn+1) have 
the same sign. (We need this fact so that we can replace the absolute value 
of their sum by the sum of their absolute values.) 


have opposite signs 


|ba pe a| = \(Tdn AF SQn41)£ = (rpn ete SPn+1)| 
= Ir (Qn — Pn) ate 8(Qn412 — Pn+1)| 
= Ir|lanz — pr| + |s||@n412 — Pn+i| 
> |r|ldnz — Pn| = |dnz — pr] 


as claimed. | 


As a consequence of this result we can now spell out the sense in which the 
convergents are best approximations to the irrational number. 


Corollary Convergents are best approximations 


Pi P2 
n° q2 7 
for any rational number 5 with 1 <b < qn 


,--- be the convergents of the irrational number x. Then, 


the equality holding only if a = pn and b = qn. 


12 
Thi f h 
is means for example, that as 165 


rational approximation - to e with 0 < b < 465. 


is a convergent of e there is no better 
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ie 
Remember 5 is in lowest terms. 


Proof of the Corollary 


a 
Assume, to the contrary, that there is a rational number 5 with 1<b< qn 


and 
a Pn 
oe dae 
Then 
ldnX — Pn| = Qn m Te wee > b|x—5| = [be -al 
contradicting the result of the theorem. E 


4.3 Rational approximations to 7 


Having seen that there are patterns for the ICFs of all quadratic irrationals 
and also for the irrational e, we would be lacking curiosity if we did not 
enquire about the ICF of the most famous of all irrationals, 7. 
Disappointingly, it transpires that there is no known pattern to the sequence 
of partial quotients of 7. We can see how the ICF begins by working with 
some good known approximation to 7. For example, take m = 3.141592654. 
With the aid of a calculator, the continued fraction algorithm yields the 
following. 


m= 3 + (m—3) = 3 + 0.141592654; a =3 
1 
eee + ey k 2 » = 
lia a2 =7 
i 
0.0625133059 ~ 15.99659454; az = 15 
mp, a is ° = 1 
099650454 ~ 1008417099; as 
1 
= — 999 6458629: es 
0.003417009 ~ -” 622; as = 29 


This, of course, is only giving a crude idea of how the ICF of m begins. Here 
we are not calculating the ICF of z itself but rather the continued fraction of 
the rational approximation to nine decimal places with which we began, and 
as we are not calculating the reciprocals exactly, rounding errors are 
creeping in. For example, if we were to use a calculator which works only to 
eight decimal places the ensuing value of a; would turn out to be 
significantly different from 292. But the partial quotients we have obtained 
thus far are correct for 7. 


* = 8,7, 15,4, 982. .sI- 
The resulting sequence of convergents begins 
3 22 333 355 103993 


The search for good rational approximations to 7 has a long history. With 
the aid of today’s computers the value of 7 is now known correct to millions 
of decimal places. The first known rational approximations can be traced 
back to Archimedes (287-212 BC). Archimedes utilized the fact that the 
circumference of the circle lies between the total edge length of any inscribed 
polygon and that of any circumscribed polygon. By using a polygon of 96 
sides he discovered lower and upper bounds given by 

223 22 

71 << 7 

22 : ; 

The value —, known as the Archimedean value of 7, is commonly used as 
an approximation to 7, being accurate to within 0.001264. The lower bound, 


Calculations like these are very 
sensitive to rounding so your 
calculator may well produce 
different results from ours. 
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2 22 
7° is not a convergent of 7 but turns out to be more accurate than —, 


333 
the error being approximately 0.0007476. The next convergent of 7, 706’ 
must be more accurate still because we know it is the best approximation of 


all rational numbers with denominator not exceeding 106; the error in 106 is 
approximately 0.000083. 


355 
However the next convergent, ——, is remarkably accurate. (This value for 7 


seems to have been first known to a Chinese mathematician Tsu Chung-Chih 
(430-501 AD) but was not discovered in Europe until late in the sixteenth 
century.) Because of the large partial quotient 292 there is a huge jump in 
the values of numerator and denominator following the convergent T13’ and 
the corollary to Theorem 4.1 gives the bound on its error: 

355 1 3 


~ TB < i97 


= 113 x 33102 ` 107° 


355 
Our theory tells us that I3 is a better approximation than any other 


rational with denominator not exceeding 113. In fact this claim is not very 
flattering, as the following example will show. 


Example 4.2 


355 
Show that a3 is the best of all rational approximations to m with 
denominator not exceeding 10000. 


Suppose that : is a rational number with 114 < s < 10000. Then 
Tie D 
U o: 
But we also have 


_ 355 
113 


_ [113r — 355s| 5 1 = 1 
is 113s =s 13510000" 


1 


$ 113 x 33 102° 


355 ., Teas 1 i 1 
Hence T13 lies within T13 x 33 102 of m but is further than 113 x 10000 


7 
from —. The triangle inequality gives 
s 


. el tis 
1 1 


> 713 10000 113 x 33102 
23 102 1 


113 x 10000 x 33 102 j 113 x 33102; 


if aI -| = 


355 
Thus - is not as good an approximation to 7 as is 13° + 


Problem 4.3 

Let x = [0, 1, 2,3, 4,5,6, ...]. Determine the first convergent of x which is 
within 1075 of z. Show that of all rational numbers ae where the positive 
denominator s does not exceed 100, the most aodai approximation 


t is. 
OF Be 
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The triangle inequality for real 
numbers says 


ja — b| < ja — z| + |x — bl. 
The form used here is 
ja — z| > |a — b| — |x — b]. 


We have seen that the convergents to the irrational number g satisfy 


peni and, as the denominator sequence {qn} is increasing, 
dn QnQn4+1 
g= à = = A question we might ask concerns whether or not this 
n 
property ‘piscina convergents; that is, if a rational number — © has the 
property |z- TE =a is 5 necessarily a convergent of x? It is not difficult 


3 1 3 
to show that the answer is no. For example, |/3 — 3 <= Z and yet 3 is not 
a convergent of /3. Nevertheless there is an inequality of this kind which 
does characterize convergents. 


Theorem 4.3 Which rationals are convergents? 


If the rational number í satisfies 


e-S]< 3 
7 ~ 5) > 20 


a, 
then A is a convergent of x. 


Proof of Theorem 4.3 
Suppose to the contrary that = is a rational number (in lowest terms) 


satisfying the given inequality but which is not one of the convergents = of 
n 


x. Let r be the unique integer for which q, < b < qr+1. Then, by 
Theorem 4.2, 


a 
lave — pel < |b —a| = 62-2] <dx 55 = a 
Therefore 
ar |e- 2 R E RS pe : 
2b’ 2q,b° 
The triangle inequality then gives 
a Pr Pr | a | 1 
bs Rema n ler ht eked: at OR : 
S alot el Obl gk (*) 
On the other hand 
Cees = qra — prb ia bes) 
b Q qrb qrb 


(The final inequality comes from the fact that qra — prb is an naa and is 
a 

non-zero, for otherwise 5 would be equal to the convergent — r which we 

have assumed not to be the case.) Combining (*) and (**) we have 


Diet a 
gb-  2qrb 202 


and this simplifies to qr > b, giving a contradiction. a 
Problem 4.4 

igg 2+vV11 
Without calculating the ICF show that TA is a convergent of aa 
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ADDITIONAL EXERCISES 


Section 1 


1 Express each of the following rational numbers as an FCF. 


(a) = (b) = (c) = (d) 3.142 


2 (a) Ifx>1 and z has FCF [a1,a2,a3,... an] show that Š has FCF 


£ 
[0, a1, @2,a3,...@,]. What is the corresponding result when 
Ori <a 


(b) Suppose that —1 < x < 0 and that x has FCF [—1, a1, a2,... an]. 


Show that [a1, a2, ... an] has value er 


3 Evaluate the sequence of convergents for each of the following FCFs 
and, in each case, verify that Theorem 1.4 holds. 


(a) [2,2,4,2,4,2]  (b) [1,2,3,4,5,6,7] 


4 Use continued fractions to find the general solution of each of the 
following linear Diophantine equations. 


(a) 1727+ 91ly =1 (b) 16x + 219y = 17 


Section 2 


1 Evaluate the following periodic continued fractions. 
(a) [(5,1,3)] (b) [0,(5,1,3)] (c) [1 2, (5,1,3) 


e2/k _ 4 r 


2 Euler showed that the ICF of the irrational number ———— i 
e2/k +1 


S 


[0, k, 3k, 5k, 7k, 9k, . ..]. 


Use this information to find values correct to four decimal places for the 


following. 
e—1 er =1 
b 
(a) = er (b) 3 =a 
Section 3 


1 Determine the ICFs of the following irrational numbers. 


(a) vil H aa (c) srg 


2 If0<k< 2nand k divides 2n, show that 


Vn? +k = fn. (Zany) 


k; 
Hence find ICFs for (a) v12, (b) v24. 


3 Determine the ICF of vn? — n, for a general n > 2. Use this to write 
down the ICFs of v2, V6 and v12. 
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Section 4 


1 


The irrational number x has ICF (0, 2, 4,6, 8, 10,...]. Find the best 
rational approximation 2 to x with: 


(a) b< 100; (b) b< 1000. 


Theorem 4.3 tells us that any rational number ; satisfying 


5 <a 


a 
is a convergent of x. It is not true that every convergent 5 of x satisfies 
the inequality. 

Pn+1 


dn dn+1 
at least one of them does satisfy this inequality. Hint: Assume the 


contrary and use the fact that 


3 p 
However, prove that of two consecutive convergents — and of z, 


Pn+1 ES Pn 
Qn+1 dn 


Pn+1 
dn+1 


Pn 
Qn 


= |T 


+- 


A theorem of Hurwitz (1891) shows that for each irrational number x 
there are infinitely many rationals ; satisfying 
p-il< 7 
z- -=| < —. 
bl ` /5b? 
(a) Why must any rational number = satisfying the above inequality 
be a convergent of x? 


(b) Find two convergents of V3 which satisfy this inequality, and two 
which do not. 


Challenge Problems 


1 


In applying the Continued Fraction Algorithm to determine the ICF of 
y/n we obtain the sequence of equations 


Vn = 2 = lan £2] = [a1, a2, 73] = [a1, a2, 03,24] =---, 


where, for all k > 1, ap = int(x,) and the step from one stage to the 
next is accomplished by 


Tk = Qk + s 
Tk+1 


Prove that for each k > 1, 


Tkt 
= 


Tk 


where {rọ} and {sx} are sequences of integers given by 
2 
a aE 
rı =0, s1=1, Tk41 =@kSk—Tk, Ski = ao ig 
k 


Use this sequence to determine the ICF of v14. 


Let x = (a1, a2, a3, a4,...]. Find the ICF of —z. 


41 


3 Let x = [a1, a2, a3, a4,...]. Show that provided az is even, the ICF of 2x 
satisfies 
2[a1, a2,43,Q4,.-. | = [2a1, tag, 2[as, Q4,45,-- a l; (*) 
and deduce that provided all the even numbered partial quotients are 
even, 
2[a1, a2,93,Q4,..- | = [2a1, tag, 2a3, tas, 285; $6, a J; 


Find a formula corresponding to (*) for the case of a> being odd. 


4 Let An = (9 : ) and M, = ( 2 ie ), where 2, 22, 
l an dn+1 Pn+1 q Q 
P3 


a ... are the convergents of [a1,a2,a3,a4,...]. Show the following. 
3 


(a) Mn = Any1Mn-1, forn>1 

(b) My = Anir AnA i. AA 

(c) By considering the determinant of M,,, deduce the formula of 
Theorem 1.3 property(a): 


Paini -aapna = (—1)"; form > 2. 


(d) By considering M,, and its transpose, deduce the formulae The transpose of ( i) 
 .. 
Daag > (ine Ont» anas 425 01] is (4 Ji 
and 
q 
= TE EE An—2,---,43, ag]. 
Qn-1 


SOLUTIONS TO THE PROBLEMS 


Solution 1.1 


1 4 30 
(a) 2,34) =24+—— =24 — == 
re BS 
+- 
4 
1 ee 
(b) [0,2,2,2] = 0 + p—-=0+— 3 =04+-5=5 
— i 2+ 
2 aa 
a 


Solution 1.2 
(a) 11=1x7+4 


E aae BS: 
4=1x3+1 
nA ec 


So = = [1,1,1,3]. 
(b) 43=2x19+5 


19=3x5+4 
§=1x4+1 
rE 


43 
So = =a LA] 
0 = 6,3,1,4 
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Solution 1.3 
(a) 65 = 2 x 23+ 19 


23=1x19+4 
19=4x4+3 
4=1x34+1 
ISSN 
65 
So TT 24a]: 
(b) 23 =0 x 53+ 23 
53 = 2 x 23+7 
23 =3x7+2 
ESS S2 
A 
23 
So 53> [0X25 35352): 
(c) —94 = —8 x 13+ 10 
13=1x10+3 
10=3x3+41 
SOO 


—94 
= pe ee 
So = [-8, 1, 3, 3] 


Solution 1.4 
(a) [2,.2;3}'= [1 2,2,4) 
(b) [—4, 3, 1,6, 1] = [—4, 3, 1, 7] 
(c) [0] = [1,1] 
Solution 1.5 
(a) The numerators of [2, 3, 4, 1,2] are 
pi=2; po=2x3+1=7; pp =4x7+2=30; 
pa=1x30+7=37; ps =2 x 37+30= 104. 
The denominators are 
= Go es = x: 
qg =1x13 +35 i6; g,=2 x 16+ 13 —45. 


The convergents are therefore 
2 7 30 37 ss 104 
L a want Gy 45° 
(b) The numerators of [—2,3, 1,6] are 
pi =—2; po=(-2)x3+1=-5; 


pa = 1x (—5)+(-2)=—7; pa =6 x (—7) + (-5) = —47. 


The denominators are 


q=]; @=3; g=1x34+1=4; qa =6x44+3=27. 


The convergents are therefore 


-2 5 -T ai 
i 2 D 


Solutions to the Problems 


Note that the integer part of -5 
is —8 (not —7). 
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Solutions to the Problems 


Solution 1.6 


(a) 


The convergents are 
43 = 18" 22 —151 
20 E E TT 


The convergents are 
Ie 2 20. Fl 26 q 71 


A r dw pa 
Solution 1.7 
(a) The convergents of (2,3, 4, 1,2] were found to be 
2 AOI ON 104 
irn mp a S 
The FCF [2,3, 4, 1,2, 2] has the same first five convergents plus a sixth, 
which is its value. The latter is determined from the previous 
convergents as 
2 x 104+37 245 
2x45+16 106° 
(b) Similarly [2, 3, 4, 1,2,100] has the same first five convergents as 
[2,3,4, 1,2] and has final convergent 
100 x 104+37 10437 
100x45+16 4516" 
The addition of a further, even larger, partial quotient has not changed 


the value of the FCF by much. In fact ~ ESTL, 


245 10437 
— =2. STL OU: 
106 3113208 and 7516 2.3111160 


Solution 1.8 
From the Euclidean algorithm 
48=2x19+10 
19=1x10+9 
10=1x9+1 
Sege 


48 
and so 19 has FCF [2,1,1,9]. Its convergents are 


From the final two convergents we have 
48x2-—5x19=1 

and therefore, multiplying through by 4 and rearranging, 
48 x 8 — 19 x 20 = 4. 


This gives x = 8, y = —20 as one solution of the linear Diophantine equation 
48x + 19y = 4, and so the general solution is 


zr=8+19t, y=—-20—48t, teZ. 


Solution 1.9 
The convergents of [2,1,1,2,1,1,2] are 
2-3. 9213. AS 80 
LP? eee H 
Theorem 1.4 orders the convergents as 
go * 
oT le aie ey j 


which is readily confirmed by evaluating these as decimals. Our results are 
recorded to a maximum of three decimal places. 


L25 < DST L2 SA 2 S. 


Solution 1.10 


For any FCF all partial quotients (with the possible exception of a1) are 
positive. So qı = 1, q2 = a2 and q3 = azaz + 1 satisfy 1 = qi < q2 < q3- 


Heading for an induction proof, suppose that qk > qk—-1 > 1 for some k > 3. 
Then 


dk+1 = Ak+14k + qk-1 = qk + Qk-1 Z qk +1> 4K, 
and the induction step is established and the proof is complete. 


Finally, since q3 > 2, the above inequality qk+1 > qk + qx—1 Shows that, 
from q4 onwards, the q,s increase in steps of at least 2 and as q4 > 3, qk = k, 
for k > 5 follows. (This could be proved slightly more formally by using 
Mathematical Induction.) 


Solution 2.1 
The convergents of the ICF [1,1,1,1,...] begin =, =, =, 5, =) => yo: 
Le ae, oe nded above b i =0:025 
ITO =r ou above by =~, = 0.025. 
1 1 
The error in Cg = z is bounded above by an 0.009615 (approx). 


The actual error in Cs is approximately 1.61803 — 1.6 = 0.01803. 
The actual error in Cg is approximately 1.61803 — 1.625 = —0.00697. 


Solutions to the Problems 


To avoid approximate division note 
that, for example, 5 x 31 < 12 x 13 


from which — < = follows. Each 


of the inequalities can be tested in 
this way by cross-multiplication. 


In fact, the ‘smallest’ sequence of 
denominators occurs when every 
partial quotient (after a1) takes its 
minimum value of 1, and then the 
sequence of denominators is 1, 1, 2, 
3, 5, 8, ..., which is the Fibonacci 
sequence. 
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Solutions to the Problems 


Solution 2.2 
If we put x = [3,1,3,3,1,3,3,1,3,...] then, as in Example 2.1, 
x = [3,1,3,a]. The convergents of [3, 1,3, 2] are 

a 4 15 1l5z +4 

oss ae 
and, as the final convergent is x itself, 

15r +4 

EATE 
As x > a, = 3, it is the positive root of this quadratic equation and 
7+ V65 
Sen” 


that is 4z? — 147 — 4 = 0. 


So 7 = 


Solution 2.3 


2 
Putting z = [2,2,2,2,...] we have x = [2,2]. The convergents of [2, x] are i 


2 1 
and sae 


a 
quadratic z? — 2x — 1 = 0 which has roots x = 1+ V2. As x > a; =2 it is 
the positive root of this quadratic equation and so 


(2,2,2,2,...]=2=14 v2. 


Now notice that we can add an integer to the value of an ICF by adding it 
to the leading partial quotient a;, that is 


. As the second of these is equal to x itself we reach the 


[a1, a2, a3,...] +n = [a1 +n, a2, 43,...]. 
Here, adding —1 to both sides, 
[1,2,2,2,..] = v2, 


which gives the ICF of v2. Its convergents begin: 
1 3 Lo Troa 99%, 22.39 
1’ 2° 5° 12° 20° 70 169° 
Now V2 lies between each pair of successive convergents. So, for example, it 


lies between 29 and 70’ whose difference in size is 30x70" Hence 39 lies 


(whic is less than l ) of V2. 


ithi 
within 1000 


1 
29 x 70 
Solution 2.4 
(a) Let x = [(1,1,1,4)] = [1,1,1,4,2]. The convergents are 
je Yeo a |! 147 4+ 3 
hata ae = 
i led Bie! 9r + 2 
The resulting quadratic equation is 9x? — 12x — 3 = 0 which, after 
cancelling the common factor 3, becomes 3x? — 4x — 1 = 0. This has the 


2 7 
positive root x = wag So 
(2541, 1)1,4))=2 a ee : 

’ ’ ? 3 ed £ F V7 

T+2V7 V7 x (V7 +2) 

24+ 77 24+ V7 

1 
(b) Let x = [(3)] = [3,2]. Then 2 = 3+ a giving the quadratic 
34+ 713 

x? — 3x — 1 = 0 which has positive root x = ee Now 


i225 a h = 42, 22] 
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2 5 5r+2 
=,= th 
has convergents 9 and FSU and therefore 
3+v13 2 
eat Maes E N 


[2, 2, (3)] 


ae oe 34/13 “$2373 
Se eee 


and this, when simplified by multiplying top and bottom by 8 — 2/13, 
gives 


pag- E 


Solution 3.1 
(a) xı = V3 lies between 1 and 2 and so int(V/3) = 1. Therefore 


V3 =1+(v3-1), giving a; = 1 and 22 = 5 — 

Repeating for x2 we have 

ag = a1 A, giving a = 1 and zs = 
From z3: 

a a 
V3-1 V3-1 


As z4 = £2 we are now in a repeating cycle with a4 = a2, a5 = a3, 
ag = a4, etc. So V3 = [1, (1, 2)]. 


(b) zı = 24+ VT = 4+ (V7 — 2); a, =4, = 
to — 1 E A =N 3 — 3 
V7-2 V7+2 3 a 3 Vil 
z3 = : VE _ az = 1 LA z 
Vi-1 vī+1 2 Ze v7-1 
%4 = 2 MER e gaea =l T5 = 3 
V7-1 Viri 3 SE i YUL 
n= x A Vima 


Hence as = a1, ag = Q2, a7 = Q3, etc. 


Therefore 2 + v7 = [(4, 1,1, 1)]. 


Solution 3.2 
For n > 1 the integer part of Vn? + 2 is n and so we write 


Vn? +2=n+(Vn?+2-n), 
iving a; = n and z2 = : 
multiplying top and bottom by Vn? + 2 + n: 


. This expression for x2 is simplified by 


1 WEFAN TAs 
n?+2—n 2 2 ; 
giving a2 = n and x3 = poe PES 
Vn? +2-n 
a n24+24+n=2n+ (Vn? +2-n), 


giving a3 = 2n and x4 = T2. 


Solutions to the Problems 


In each of the last two calculations 
we multiplied top and bottom by 
V3 + 1 and then pulled out the 
integer part. 


47 


Solutions to the Problems 


This pair of equations now cycle and so Vn? + 2 = [n, (n, 2n)]. 


Hence we have: 


v3 = (1, (1, 2)]; 
V6 = [2, (2, 4)]; 
Vi = [3, (3, 6)]; 
V18 = (4, (4,8)]. 


Solution 3.3 


(a) From the evidence presented in Table 3.1 you might have made some of 
the following observations: 


e the ICF of yn is periodic, with cycle starting at av; 
e the final partial quotient in the cycle is 2a;; 


e leaving aside the final term 2a,, what remains of the cycle is 
symmetric (palindromic), reading the same forwards and backwards; 


e the palindromic part sometimes has even length, sometimes odd, 
and in some cases (e.g v2 and v5) is non-existent. 


(b) The effect is as follows. 


Adding a, which is equal to the integer part of y/n, gives a first partial 
quotient of 2a;, and this can be absorbed into the cycle leaving a purely 
periodic ICF. For example, 


VT + int(V7) = [(4, 1,1, 1)]. 


(c) Collecting all this together we might make the conjecture that the ICF 
of \/n takes the form 


vn = [a1, (a2, a3, sv . a3, a2, 2a1)]. 


Solution 4.1 
The convergents of v5 are 


2 9 38 161 
ete T eS 


The convergent = satisfies 


2x17x72 |2 17| <I7x 72° 


That is, working to five decimal places 


1 z| 1 


0.00041 < 


38 
Ee 0.00082. 
=e. 


Solution 4.2 
The convergents of v2 are 


13 7 17 41 99 239 
2" 6°43" 2a" Wp? des" °° * 


99 
The convergent 70 satisfies 
99 1 1 
ees 


70| Č 70x 169 ~ 10000 


and so has the required accuracy. It is the first convergent with this 


= 0.0001, 


property because the previous convergent 29 satisfies 


z = 0.0001. 


41 1 3 1 
29 2x 29x 70° 10000 
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Solution 4.3 

The convergents of x are 
w t-22- T7 WI R2 
Toe s a0 oe" A 


Pe a 
The convergent 725 satisfies 


T 


157 1 
225 225 x 1393 


and so has the required accuracy. It is the first convergent with this 


< 0.000004, 


: 30 z 
property because the previous one, B’ satisfies 
30 
ae 3 aaa AAA 


Let be a rational different from 5 with 1 < s < 100. Then, since 


43r — 30s is a non-zero integer, we have 


m - ee 1 
s 43 43s = 438 43x100 
But 
ae SO ee 
43 43 x 225’ 
and so the triangle inequality gives 
y: i 1 125 
-7| > 13x100 43x225 43x 100 x 225 
i 30 
? 7B x 225 ara 


(ON sai r 
and so 73 88 better approximation than = 


Solution 4.4 


17 
By Theorem 4.3, a sufficient condition for E to be a convergent of 
2+v11 17 it 
is that Ere : 
5 16) 3216 


This can be confirmed by use of a calculator or 


2+v11 17|  |—53+16y11| 


5 p 80 


|(—53 + 16V11)(—53 — 16V11)| 
80 x |(—53 — 16v11)| 
7 7 
-~ 80 x (53 + 16/11) 
a =e < : $ 
8000 ~ 2 x 162 


Solutions to the Problems 


2+vil 
5 


537 = 2809 
11 x 167 = 2816 


53+ 16/11 > 53 + 16 x 3 > 100 
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SOLUTIONS TO ADDITIONAL 
EXERCISES 


Section 1 


1 (a) 23=0x 54423 


54=2 x 23+8 
23 =2x8+7 
S—1oGy ol 
TPT 


23 
So 54 S 2,2, 1, 7). 
(b) 117=6 = 17 +15 
i A a a 2 
LS E A R 
2 De 
117 
So 7 = A]: 
(c) —37 = (—2) x 24411 
24:2 xa ED 
L596 A 
i eel Gu | 


Geg [=2;9.5; 9). 


24 
3142 
.142 as ——. 
We treat 3 as 1000 


3142 = 3 x 1000 + 142 

1000 = 7 x 142+6 

142 = 23 x6+4 
6=1x4+2 
eae AS Ga, 


So 3.142 = [3, 7, 23, 1, 2]. 


(ad 


wa 


2 (a) If x = [a1, a2, a3,...an] > 1 then 
1 1 1 
—~ =04+—=04+——_———. = (0.41, a9, Gg, Ga: 
£ T [a1, a2, a3, .. . an] 
Note that a; is not the leading partial quotient in this FCF and so 
it is essential that a; > 1, a fact which is guaranteed by zx > 1. 


1 
For 0 < x < 1 we have > 1, and so we can apply the result just 
1 1 
obtained to —. If =e [b1, b2, b3,- . bn] then the FCF of its 
x 


reciprocal (noting that E x) is given by x = [0, by, bo, b3,... bn]. 

In other words, if 0 < x < 1 then its FCF will have leading partial 
1 

quotient 0, say x = [0, a1, a2, . . . an] and then A [a1,a2,...An]. 


(b) If x = [—1, a1, a2, ... an] then 1 + 2 = [0,a1,a2,...an]. By the 
result of part (a) the reciprocal of 1 + x has FCF [ay, ag,... an]. 
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Notice that this works without 
reducing the rational numbers to 
lowest terms. We leave you to 
ponder why. Generally we consider 
our rational numbers to be in 
lowest terms. 


3 (a) [2,2,4,2,4, 2] has convergents 
2 5 22 49 718 485 
T 2 9° 20° 89 i 198 
According to Theorem 1.4 we have the inequalities 
1 9 89 108-20. 2” 
each of which is readily verified. 
(b) [1,2,3, 4,5, 6,7] has convergents 
1.3 W G 22 13 9976 
Kr roe Goel 
The inequalities given by Theorem 1.4 are 
EE OE E NR a 
T «0: G06) ~ 872 ~ 30 ~ 2’ 
each of which is readily verified. 


91 
4 (a) The FCF of 7 is [5,2,1,5] and it has convergents 


p- ai 46 91 
Ton Sena i 
From the final pair of convergents 3 x 91 — 16 x 17 = 1, which 


shows that + = —16, y = 3 is one solution of 17x + 9ly = 1. The 
general solution is therefore given by 


x=-16+91t, y=3-17t, tez. 


219 
(b) The FCF of Je is [13, 1,2,5] and has convergents 


~ 


13 t4 41 219 

7° ee and 6" 
From the final two convergents we have 3 x 219 — 16 x 41 = 1. 
Multiplying through by 17 gives x = (—41) x 17 = —697 and 
y = 3 x 17 = 51 as one solution of 16x + 219y = 17, and the general 
solution is given by 


x = —697 + 219%t, y=51-16t, tez. 
Alternatively since t = 3 gives the solution x = —40, y = 3, 


«=—40+219, y=3-16t, TEZ 


Section 2 


1 (a) Let x = [(5,1,3)]. Then x = [5,1,3,2] and this has convergents 
5 6. 23 23z + 6 
=, >, — and ———. 
seat ps! AF 
As the final convergent is equal to x we have the quadratic equation 
4x” — 227 — 6 = 0. 


11+ v145 


The positive root of this equation is x = 3 
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(b) 


m, 
i) 
~ 


1 
(0, (5, 1, 3)] has value = where z is as in part (a). That is 


4 
0, 5,1,3 =- a) 
[0,4 2 114+ v145 


and this can be simplified by multiplying top and bottom by 


11 — V/145: 
[0, (5, 1, 3)] = = z dei 


Is 
2 E= -,- : 
[1, 2, (5, 1,3)] = [1,2, £] and has convergents > and ap Pl 
substituting for x from (a), 


_ 37437145 
26 + 2/145’ 


23 + V145 
24 : 


[1, 2, (5, 1, 3)] 
which can be simplified to 


cid 
2 = (0,2,6, 10, 14, 18,...]. The 
e+1 


Taking k = 2 we have 
convergents begin 
0 1 6 61 860 
Ya 1a eee 8 


61 1 
Th . ee Beige 
e convergent 132 $ accurate to within 152 x 1861 which is less 


than 5 x 1076 and guarantees correctness to four decimal places. 


2 
Taking k = 1 we have s 7 ; = [0,1,3,5,7,9,...]. The convergents 
are 
O T3 rt6 T5- 0a: 
Vt Soa es oe 
The convergent ne is accurate to within ee which is less 
151 151 x 1380 


than 5 x 10~® and guarantees correctness to four decimal places. 


Section 3 
1 
1 a) 2 =v11=3+(v11—-3); aj = 3, T= 
(a) sı =v =3+(Vi1—3) 1 Bs roa 
1 VIFS fit =2 2 
= =e ee = 
viL— 3 2 2 viļ-3 
a = V114+3=6+(V11-3); a3=6, 2 =——— 
As £4 = X2 we are now in a cycle with the second and third 
equations recurring indefinitely, so v11 = [3, (3, 6)]. 
1+ 713 V13 —3 2 
b ry = =2+ ; a =p x= ee 
(b) 2 J 7 1 2 VERE 
2 V13 +3 vis—3 2 
= =3 = i an = 
JB=3 2 2 V13—3 


52 


As £3 = £2 we are immediately into a cycle and, after the initial 
partial quotient of 2, the partial quotient 3 recurs indefinitely. 


1+ 713 


Hence n [2, (3)]. 


2 


(c) 5 SN Met aa ee 
1 5 = 5 ’ 1 , A ae 
a Fie ee E E 
a, a oe 4 EE 4 ’ 2 ’ ar at ay 
BOS a en oy aaa 
Jn=3 . 3 = 3 , 3 = 4; S 
> yarr ey oe 
aA 4 4 ’ 4—4, 5 Jit! 
$e oT ee See 
vaL—1 5 5 ’ 5 , 6 v2 —4 
= /214+4=84(V21-—4); ag =8, oo 
eg VA ); a6 = 
ee et ge 
Jii— 4. 5 = 5 ? 7 ? anes 


Now zg = z2 and we are in a cycle, with the sequence 1,2,1,1,8,1 
of partial quotients recurring. So 


4 21 
= ZE AT ieee a 
But notice that this gives 
2 
2 = = [(1,1,2,1,1,8)]. 


In fact, although it was not immediately obvious, £7 = 2; and we 
were in a cycle at the point of calculating x7. 


Since 0 < k < 2n we have n? < n? +k <n?+4+2n+1=(n+1)?. 
Therefore the integer part of Vn? + k is n, and so 


1 
Vn? +k=n+(Vn?+k-n); a =n, t2 = ————, 


n2? +k-n 
1 _ VRE Rn 
Vn? +k-n k ` 
The integer part of Vn? + k +n is 2n, and as k divides 2n we have 
1 Qn Vn?+k—-n 2n k 
= ate ; 2=—, 23 = -= > 
m+k—-n k k k ne+k—n 


k 
IY ol PRT eel = ia a Da as = 2n; 


1 


Pioa eee 
ia Vn? +k-n 


2 
As z4 = z2 we are in a cycle with the pair of partial quotients oa 2n 
repeating indefinitely. So 


Vn? +k =- In (F, 2n)] } 


k 


(a) Take n = 3 and k = 3. Then V12 = [3, (2, 6)]. 
(b) Take n = 4 and k = 8. Then V24 = [4, (1,8)]. 
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3 For n> 2 we have 
(n—1)? =n? -2n4+1<n?-n<n? 
and so Vn? — n lies strictly between n — 1 and n. Therefore 
w= Vn?—-n=n 1+ (Vn? —n-(n—1)), giving a =n—1. 
Continuing the Continued Fraction Algorithm: 


g- NNO Taen a 
ana na a 


As the integer part of Vn? — n is n — 1 the integer part of 
Vn? -n+n-1 


is 2 and so 
n— 1 
Vn? —-n—=(n-1 
z2 =2+ é = - ) giving a2 = 2. 


At the next stage we look at 


asi 
s= ey V e+e 


As this has integer part 2(n — 1), 


z3 = 2(n — 1) + (Vn? -n — (n — 1)), giving az = 2(n — 1). 


And there is good news. The fractional part, Vn? — n — (n — 1), is the 
same as the fractional part of 21, so we are in a cycle and can conclude 
that 


Vn? -=n = [n—1, (2, 2n — 2)). 


For example: 
n = 2 gives V2 = [1, (2, 2)] = [1, (2)]; 
n = 3 gives V6 = [2, (2, 4)]; 
n =A gives V12 = [3, (2, 6)]. 


Section 4 


1 (a) The convergents of x = [0,2,4,6,8,10,...] are 
0 1 4 25 204 2065 
12-97 NOP RR? ang Abo: oe 

1 


2 
-56| ` 56 x 457 
be the best rational approximation amongst rationals with 


T and is known to 


25 
The convergent 56 satisfies 


denominator not exceeding 56. Now suppose that $ is any rational 


2 
with b < 100 and different from =. Then 


ee teori 


b 56 56b © 56b ^ 56 x 100° 
It follows (from the triangle inequality) that 
a ie 29 25 
i-42-8]--8 
b b 656 56 
1 1 357 1 


> 56x100 56x457 56x 100x457 56x 457’ 


25 
and therefore 56 is a better approximation than = 


54 


204 
(b) In the same vein, the convergent 157 satisfies 
204 1 
*~ 457 | © 457 x 4626" 


4 
So if ; is a rational number different from = and with b < 1000 


457 
then 
eee 1 
b 457 457b = 457b ^ 457 x 1000’ 
and 
a 1 1 3626 1 
P S z| > 757 x 1000 457x4626 457 x 4626 x 1000 457 x 4626" 


204 
Therefore, —— is the best approximation amongst rationals with 


denominator not exceeding 1000. 


Since the convergents of x are alternately greater than and smaller 


than z, 
Me Pn+1 m fe _ Pr x Pnt1 — Pn = \Pn+19n — PnQn+1| x il 
dn+1 dn dn+1 dn QnQn+1 GnQn+1 
So if neither convergent satisfies the given property 
1 il 1 


QnQn4+1 : 202 44 i 242 
This can be rearranged as (qn+1 — qn)? < 0. As a square cannot be 
negative and qn+1 > qn except for the possibility of qı = q2 = 1, this 
gives the required contradiction for all except this exceptional case. 
However, when q1 = q2 = 1 the first two convergents of x will be int(x) 
and int(x) + 1 and the result makes the modest, and true, claim that 


one of these lies within 3 of the true value. 


(a) Since V5 > 2, |e G and Theorem 4.3 guarantees 


ERF 

bl ` \/5b2 ` 26?’ 
a, 

that z is a convergent of x. 


(b) v3 = [1, (1, 2)] and has convergents 


125 7 19 26 7 
I 1 3e a ie E 
1 1 
It can be checked by calculator that |\/3 — 1| > —, |V3—-2| < —, 
5 1 ti 1 1 5 
3— -| > —- and |v3-— -| < ——-, so the convergents — and +- 
vs 3|° 9v5 uF 4| ` 16/5 z 1 3 


7 
fail to satisfy the inequality while the convergents I and 1 do satisfy it. 
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